Environment-assisted tunneling as an origin of the Dynes density of states 



O 

(N 

< 

o 
o 

u 



i 

C 

o 
o 



> 
m 
in 
oo 
cn 

o 
o 



X 



J. P. Pekola, 1 V. F. Maisi, 2 S. Kafanov, 1 N. Chekurov, 3 A. Kemppinen, 2 

Yu. A. Pashkin, 4 0.-P. Saira, 1 M. Mottonen, 1 ' 5 and J. S. Tsai 4 

'Low Temperature Laboratory, Aalto University School of Science and Technology, P.O. Box 13500, FI-00076 AALTO, Finland 
2 Centre for Metrology and Accreditation (MIKES), P.O. Box 9, 02151 Espoo, Finland 
'Department of Micro and Nanosciences, Aalto University School of 
Science and Technology, P.O. Box 13500, FI-00076 AALTO, Finland 
* NEC Nano Electronics Research Laboratories and RIKEN Advanced Science Institute, 
34 Miyukigaoka, Tsukuba, Ibaraki 305-8501, Japan 
5 Department of Applied Physics/ COMP, Aalto University School of 
Science and Technology, P.O. Box 15100, FI-00076 AALTO, Finland 

We show that the effect of a high-temperature environment in current transport through a normal 
metal-insulator-superconductor tunnel junction can be described by an effective density of states 
(DOS) in the superconductor. In the limit of a resistive low-ohmic environment, this DOS reduces 
into the well-known Dynes form. Our theoretical result is supported by experiments in engineered 
environments. We apply our findings to improve the performance of a single-electron turnstile, a 
potential candidate for a metrological current source. 



Introduction — The density of states (DOS) of the car- 
riers governs the transport rates in a mesoscopic con- 
ductor [![, e.g., in a tunnel junction. Understanding the 
current transport in a junction in detail is of fundamen- 
tal interest, but it plays a central role also in practical 
applications, for instance in the performance of supercon- 
ducting qubits [2J , of electronic coolers and thermometers 
3, and of a single-electron turnstile to be discussed in 
this Letter When one or both of the contacts of a 
junction are superconducting, the one-electron rates at 
small energy bias should vanish at low temperatures be- 
cause of the gap in the Bardeen-Cooper-Schrieffer (BCS) 
DOS [5|. Yet, a small linear in voltage leakage current 
persists in the experiments that can often be at- 

tributed to the Dynes DOS, a BCS-like expression with 
life-time broadening XjJ, |l2j. A junction between two 



leads admits carriers to pass at a rate that depends on 
the DOS of the conductors, the occupation of the en- 
ergy levels, and the number of conduction channels in 
the junction 13[. In general, basic one-electron tunnel- 
ing coexists with many-electron tunneling, for instance 
co-tunneling in multijunction svstems |14j . or Andreev 
reflection in superconductors \wL [l6| . However, when 
the junction is made sufficiently opaque, a common sit- 
uation in practice, only one-electron tunneling governed 
by the Fermi golden rule should persist. We demonstrate 
experimentally that the sub-gap current in a high-quality 
opaque tunnel junction between a normal metal and a su- 
perconductor can be ascribed to photon assisted tunnel- 
ing. We show theoretically that this leads exactly to the 
Dynes DOS with an inverse life-time of e 2 kBT CDV R/h 2 , 
where T env and R are the temperature and effective re- 
sistance of the environment. 

We employ a tunnel junction with a normal metal- 
insulator-superconductor (NIS) structure, see Fig. QJa). 
The essentially constant DOS in the normal metal ren- 
ders the NIS junction an ideal probe for the supercon- 
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FIG. 1: (color online) (a) Geometry of the measured single 
NIS junctions made of aluminium (low contrast) as the super- 
conductor and copper (high contrast) as the normal metal. 
The tapered ends lead to 250 x 250 )im 2 pads, (b) Typi- 
cal IV characteristics, measured at 50 mK for a sample with 
Rt = 30 kQ,. Linear leakage is observed deep in the gap re- 
gion \eV\ < A ~ 200 /ieV, consistent with the Dynes model 
using 7 = 1.8 x 10~ 4 , shown by the cyan line. 



ductor DOS. Due to the BCS energy gap in an NIS sys- 
tem, the tunneling current is expected to be exponen- 
tially suppressed with decreasing temperature. Yet in the 
experiments a small sub-gap current persists as shown in 
Fig. UTtO. This leakage is typically attributed to Andreev 
current [17H20j |. smeared DOS of the superconductor [21 
non-vanishing DOS in the insulator within the gap 
non-equilibrium quasiparticles 22j , or physical imperfec- 
tions in the junction. Our junctions, like the one in Fig.[T] 
are made opaque with large normal-state resistance Rt 
to efficiently suppress the Andreev current. A convenient 
way to account for the smearing of the IV characteris- 
tics is to use the so-called Dynes model llj, [l2( based on 
an expression of the BCS DOS with life-time broaden- 
ing. The Dynes DOS, normalized by the corresponding 
normal-state DOS, is given by 
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where A is the BCS energy gap. A non-vanishing 7 intro- 
duces effectively states within the gap region, \E\ < A, 
as opposed to the ideal BCS DOS obtained with 7 = 
resulting in vanishing DOS within the gap. This model 
reproduces the features observed in our measurements 
as is shown in Fig. Hfb). We show that, effectively, 
the Dynes DOS can be produced from the ideal BCS 
DOS by weak dissipative environment at temperature 
Tem ^ A/fcs promoting photon assisted tunneling. A 
similar environment model with comparable parameter 
values has also been introduced by other authors to ex- 
plain, e.g., observations of excess errors in normal-state 
electron pumps 23, HH and Andreev reflection domi- 
nated charge transport at low bias voltages in NISIN 
structures [251 ] . 

Theoretical results — For inelastic one-electron tunnel- 
ing, the rates in forward (+) and backward (— ) directions 
through an NIS junction can be written as 
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FIG. 2: (color online) The probability density P(E) calcu- 
lated for a = 10~ 3 = -RTenvfcs /(RqA) and for a few values 
of the environment temperature, T cnv = 0.042, 0.42, and 4.2 
K. The dotted line is the corresponding Lorentzian limit of 
Eq. (j4]). The left inset shows the employed circuit model 
of an NIS junction (the rectangular symbol at the bottom 
right) in an RC environment. The right inset shows the calcu- 
lated zero-bias conductance of the NIS junction as a function 
of the capacitance C with the environment corresponding to 



f N (E T eV)[l - f S (E')]P(E - E'), (2) a = 10~ 3 and T env = 0.42,4.2, and 42 K. We use A = 200 



at bias voltage V. Here, P(E) refers to the probability 
density for the electron to emit energy E to the environ- 
ment 1 1 31 ] . The occupations in the normal and supercon- 
ducting leads are given by the Fermi functions f^/g(E) = 
^ e E/(k B T N/s ) _|_ respectively. In an ideally voltage- 

biased junction, P(E) = 6(E). The current through the 
junction at low temperature of the leads, Tn,T$ — > 0, is 
then I°(V) = e(T° + - T°_) = ^ f° V dEn s (E). Thus we 
obtain the well-known expression for the conductance of 
the junction as 



/ieV ~ ks x 2.3 K for aluminum. 



energies E <C K(RC) 1 , /csT onv as a Lorentzian [26 



P(E) 



ttA a 2 + (E/A) 2 
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d!°/dV = R^ 1 n s (eV). 



(3) 



In view of Eq. ([lj, the non-zero linear conductance at 
low bias voltages, typically observed in experiments as 
shown in Fig.QJb), suggests that the superconductor has 
non-vanishing constant density of states within the gap. 
Within the Dynes model of Eq. ([T]), the normalized DOS 
at low energies, \E\ <C A, equals 7, the ratio of the con- 
ductance at zero bias and that at large bias voltages. 
Although this approach is correct mathematically, it is 
hard to justify the presence of subgap states physically. 

Here, we base our analysis on the pure BCS DOS 
(7 = 0) and show that the Dynes model in Eq. (J]) is con- 
sistent with weakly dissipative environment when Eq. (j5J) 
is used to obtain the current I(V) = e(T + — The 
effective resistance value R — aR cnv of the environment 
arises generally from a possibly larger real part of the 
environment impedance, R cnv , which is suppressed by a 
factor of a due to low-temperature filtering, see Fig. [3] 
With this environment, one obtains the probability den- 
sity P(E) in the limit of small R <C Rq = h/e 2 and for 



where a — RkBT cnv / (RqA). As the current of an NIS 
junction is determined by the values of P(E) at \E\ < A, 
we can apply Eq. (TJ} when fc_BT cnv > A, see Fig. [2] for a 
numerical demonstration. For a general symmetric P(E) 
and Tn,Ts —> 0, one obtains from Eq. ^ in analogy 
with Eq. ©: I(V) = ^ / eV dEn°(E), and 

dl/dV = R T 1 n a s (eV), (5) 



where the effective DOS is given by the convolution 



dE'n s (E')P(E - E') 



(6) 



For the weak resistive environment described by Eq. ^ , 
the convolution of a Lorentzian gives 



rf s (E) = 5Re( 



E/A- 



(7) 



This expression is identical to the Dynes DOS in Eq. (TTJ) 
by setting a — 7, with the equivalent inverse lifetime 
e 2 kBT env R/h 2 . The correspondence between the P(E) 
theory and the Dynes model, our main theoretical result, 
is valid for non-zero lead temperatures as well, as we show 
in the supplementary material (2(|. Below, we present 
numerical and experimental studies verifying our claim. 
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FIG. 3: (color online) Measured IV curves of an NIS junction 
with Rt = 761 kQ, on the ground plane (solid symbols), and 
of a similar junction with Rt = 627 kfi without the ground 
plane (open symbols). Solid lines present the results of the 
full P(E) theory for capacitance C = 10 pF (red line) and 
C — 0.3 pF (green line). The resistance and the temperature 
of the environment are set to 7? = 2 Q, and T cnv = 4.2 K, re- 
spectively, and A = 200 /ieV. The dashed lines correspond to 
the Dynes model with the parameters yielding the best fit to 
the data. The normalized zero bias slope is 5.3 x 10~ 4 for the 
green line, and 2.6 x 10 for the red line. The inset shows IV 
curves based on the full P(E) calculation as functions of the 
shunt capacitance C. The red and green lines are reproduced 
on this graph from the main figure. 



Experiments on single junctions — The leakage induced 
by the electromagnetic environment can be decreased 
by efficient rf filtering of the leads and electromagnetic 
shielding of the sample. One way to do this without af- 
fecting the properties of the junction itself is to increase 
the capacitance C across it, see Fig. [2] In this way one 
approaches the case of an ideally voltage biased junction. 
In Fig. [3J we present the zero-bias conductance of an NIS 
junction as a function of the shunting capacitance, C, 
based on the full numerical P(E) calculation. For low C, 
the result using Eq. (J7) is valid, but for sufficiently high 
C, i.e., for ARC/h > 1, the leakage decreases signifi- 
cantly demonstrating that capacitive shunting is helpful 
in suppressing the photon assisted tunneling. 

To probe the effect of the capacitive shunting in our 
experiments, we introduced a ground plane under the 
junctions. The junctions were made on top of an oxidized 
silicon wafer, where first a conductive 100 nm thick Al 
layer working as the ground plane was sputtered. On top 
of this, a 400 nm thick insulating high-quality AI2O3 film 
was formed by atomic layer deposition. The junctions 
were patterned by conventional soft-mask electron beam 
lithography on top of AI2O3. For comparison, junctions 
were made both with and without the ground plane. 

The experiments reported here were performed in a 
3 He- 4 He dilution refrigerator with a base temperature of 
about 50 mK. All the leads were filtered using 1.5 m of 
Thermocoax cable between the 1 K stage and the sample 



stage at the base temperature. The IV curves such as 
the one in Fig. [ljb) are thermally smeared at elevated 
temperatures, but below 200 mK we observe hardly any 
temperature dependence. Figure [3] shows the IV curves 
measured at the base temperature for one junction on 
top of a ground plane and for a similar junction without 
the ground plane, together with numerical results from 
the P(E) theory. The capacitive shunting decreases the 
zero-bias conductance significantly. The shunt capaci- 
tance values employed in the P(E) theory, 10 pF and 
0.3 pF, respectively, match well with the estimates for 
the experimental values in each case. The sample with- 
out a ground plane with C — 0.3 pF is already entering 
the regime, where the capacitance is too small to play 
a role. We used an effective environment resistance of 
R = 2 f2 at T env = 4.2 K, close to the values inferred by 
Hergenrother et al. |25j in the case of incomplete shield- 
ing. However, the choice of T env is somewhat arbitrary 
here: T cnv > 4.2 K with correspondingly lower R would 
yield a slightly improved fit to the data, but T env = 4.2 
K, the temperature of the outer shield, was chosen as 
a natural surrounding in the measurement set-up. Our 
results with capacitive shunting, on the other hand, cor- 
respond to much improved shielding in the lan guag e of 



25 are 



Ref. [25[ . Although the experiments of Ref. 
quite different from ours, their situation resembles ours 
in the sense that photons with a very high frequency of 
A/ft, > 50 GHz are responsible for tunneling. 

SINIS turnstile — As a practical application, we discuss 
the SINIS turnstile which is a hybrid single-electron tran- 
sistor (SET) and a strong candidate for realizing the unit 
ampere in quantum metrology 0, [27l - [30j | . In the previous 
experimental studies ji , 27 , 2i| , its accuracy was limited 
by the sub-gap leakage. Here we test the influence of the 
ground plane on the flatness of the current plateaus at 
multiples of e/, where / is the operating frequency. The 
ground plane had a 20 /im wide gap under the SET to 
reduce the stray capacitance to the rf gate. The ground 
plane layer was covered by a 300 nm thick insulating layer 
of spin-on glass, on top of which the rf gate and dc leads 
were evaporated. Another 300 nm spin-on glass layer was 
used to cover the rf gate, and the SET was fabricated on 
top of this layer. The device is shown in Fig. Sta). This 
sample geometry is designed for parallel pumping [30j j . 
but here we concentrate on a single device. 

Figure @|b) shows that in this case, the introduction of 
the ground plane reduces the sub-gap leakage by roughly 
two orders of magnitude as opposed to a typical turnstile 
without the ground plane (the latter data from Ref. (29j|). 
In the turnstile operation, the current was recorded as a 
function of the amplitude of the sinusoidal rf drive, A g , at 
several bias voltages. In Fig.[4tc), we show the quantized 
current plateau at / = 10 MHz, and the averaged current 
on this plateau is given in Fig. 0Jd) as a function of the 
bias voltage. The differential conductance at the plateau 
divided by the asymptotic conductance of the SET is 
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2 x 10 -6 . This result is much improved over those of the 
earlier measurements [i , 27 , 28 1 , and that of the reference 
sample without the ground plane. 
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FIG. 4: (color online) (a) Scanning electron micrographs of 
the SINIS turnstiles, (b) Sub-gap IV curves of the measured 
transistors in the gate-open state (charge degeneracy). The 
slope of the linear fit corresponds to the leakage of 0.5 x 10 -6 
for the sample with the ground plane (filled circles) and 
50 x 10~ 6 for the sample without the ground plane (open 
diamonds) in units of the asymptotic conductance of each 
SET. (c) Current through the turnstile on the ground plane 
as a function of the amplitude of the applied sinusoidal gate 
drive at / = 10 MHz. The gate offset was set to the charge 
degeneracy point, and the bias voltage was varied uniformly 
between VJ, = 140 . . . 290 /iV. (d) Current at the first plateau 
as a function of Vb obtained from data similar to those in (c) 
(filled circles) showing leakage of 2 x 10 -6 and for the sample 
without the ground plane (open diamonds) with leakage of 
100 x 10 -6 and a reduced step width. 

In conclusion, we have shown analytically that the 
Dynes density of states can originate from the influence 
of the electromagnetic environment of a tunnel junction, 
and it is not necessarily a property of the superconduc- 
tor itself. Our experiments support this interpretation: 
we were able to reduce the leakage of an NIS junction 
by an order of magnitude by local capacitive filtering. 
We stress that capacitive shunting does not necessarily 
suppress the subgap leakage of an NIS junction, if the 
leakage is caused by the poor quality of the junction, or 
by true states within the gap due to, e.g., the inverse 
proximity effect Q. Protecting the junctions against 
photon assisted tunneling improves the performance of, 
e.g., single-electron pumps. Contrary to the resistive en- 
vironment aiming at the same purpose [281 ] . capacitive 
shunting does not limit the tunneling rates. 
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TUNNELING RATES IN A RESISTIVE ENVIRONMENT 

Let us consider an NIS tunnel junction with a normal-metal and a superconducing lead separated by a thin oxide 
layer, shown in Fig. 1 of the original article. We assume that the junction is in an ideal voltage bias V and neglect 
the charging energy of the single junction justified by the large capacitance provided by the leads. As in the fully 
normal-metal case sequential tunneling rates for single-electron events along/against the bias are obtained 

to be 



r± = 



i 

e 2 R T 



/OO pOO 
/ dEdE' n s (E')f N (E T eV)[l - fg(E')]P(E - E'), 
-oo J — OO 



(S-l) 



where Rt is the asymptotic resistance of the junction, ng is the Bardeen-Cooper-Schrieffer (BCS) density of states 
(DOS) in the superconductor, and the Fermi functions for the superconductor fg and the normal metal /jy are governed 
by the temperatures Tg and Tjy, respectively. The function P(E) can be interpreted as the probability density to 
emit energy E to the environment. The average electric current through the junction is given by I = e(r + — T_). 
Thus the current -voltage (IV) response of the system, which is typically the only experimental observable for a single 
junction, is given by 



i />oo />oo 

I(V) = — / dEdE'n s {E')[l- f s (E')][f N {E-eV)- f N (E + eV)]P{E-E'). 

C--^T J —oo J —oo 



(S-2) 



We note that in an ideal case with no coupling to an external electro-magnetic environment in the tunneling process, 
P(E) is a delta function and Eqs. (|S-1HS-2|) reduce into the typical Fermi golden rule results with one-dimensional 
integrals. Generally, P(E) can be written as 



cxp ( J(t) + -Et ) dt, 



(S-3) 



and by modeling the environment by an impedance Z(w) in a thermal equilibrium at temperature T env , one obtains [l[ 



J{t) 



dcu Re[Z(u)] 
uj 2ttRq 



coth 



2fcsT c: 



(coswi — 1) — zsinwt 



(S-4) 



where Rq = fi/e 2 is the resistance quantum. For a purely resistive environment, Re[Z(w)] = R/[l + (ujRC) 2 }, where R 
is the resistance of the environment and C is the total capacitance at the junction including the junction capacitance 
and parallel shunt capacitors. Although this rather simple model has essentially only two free parameters, R and T^rw, 
for low values of capacitance C, it has been successfully applied to explain several experimental observations @, S|- 
In our case of a resistive environment, the integral in Eq. (|S-4I) can be evaluated to yield 



J(t) 



cot(J3)(l 



b-h 



ill 
B 



2 £-\ 

n=l T17T 1 



-nir\T\/B 



(nir/BY 



i sign (r) (1 — e ' r ') 



(S-5) 



where p = R/Rq, t = t/(RC) and B = h/(2k B T env RC). For \t\ > 1,B, we can neglect all the terms in Eq. (jS-5|) 
containing |r| in the exponent. Use of this approximation for evaluating P(E) according to Eq. (jS-3[) is justified if the 
energy exchange E with the environment satisfies \E\ <C h/(RC), kBT eDv . For low bias voltages in an NIS junction, 
an electron needs to absorb energy of A — eV in order to tunnel. Therefore, the BCS gap gives the relevant scale of 
the energy exchange at low bias voltages. At higher voltages, the current is essentially linear in V and insensitive to 



2 



the functional from of P(E). Hence the conditions above can be taken to A <C h/(RC), kBT env . Thus, under these 
circumstances we simplify Eq. (jS-51) into 



a + — + i sign (r) 



(S-6) 



where a is the sum of all terms in Eq. (jS-31) independent of r. By evaluating the integral in Eq. (jS-31) with J(t) from 
Eq. (|S-6|) . we obtain 

p(s) = 1 c -pa/2 (pk B T cnv ) cos Qp/2) + £sin (p/2) 

{pk B T em f + £ 2 

Here, the normalization condition P(E)dE — 1 yields e _pQ / 2 = 1/ cos(/?/2). Furthermore, the situation p <C 1 is 
met in experiments @, I3| unless the environment in the immediate vicinity of the junction is deliberately engineered 
to be high Ohmic Together with the assumption \E\ <C fcs7e nv cast above, Eq. (|S-7[) assumes the form 



(pk B T ej 
* (pk B T CQV ) 2 + E 2 



p(E) = ± - , y Z B T - ■ (s- 



We note that since we neglected the term Esm(p/2) in Eq. ([5^7]) . the detailed balance P(-E) = e - E / ( - kBTelxv ) p(E) 
is not perfectly satisfied. Due to our assumption of high temperature \E\ <§; ksT env , the error due to the slight 
breaking of the detailed balance is, however, small. As shown below, the fact that the P(E) function is symmetric, as 
in Eq. (IS-81) , is crucial in reducing IV characteristics given by Eq. (|S-2p into the typical form with a one-dimensional 
integral in the energy domain. In this form, the effect of the environment is only to effectively modify the BCS DOS. 

EFFECTIVE DENSITY OF STATES 



Let us come back to Eq. (jS-21) yielding the current-voltage characteristics of an NIS junction. Since the BCS DOS 
vanishes inside the energy gap \E\ < A, Eq. (|S-2[) can be expressed in the form 

-| />oc />oo 

I(V) = — dE dE'n s (E')[f N (E-eV)-f N (E + eV)}P(E~E'), (S-9) 
eUr J-oo Jo 

if the temperature of the superconductor satisfies T$ <C A//cg. This condition is met in the experiments reported 
here since T s < 100 mK <A/fc B w2K. We expand the integral in Eq. (|S-9|) as 



o 



I(V) = — J dE I dE' n s {E')[f N (E - eV) - f N (E + eV)]P(E - E') + 



)• (s-io) 



i />OC poo 

— dE dE'n s (E')[f N (E-eV)-f N (E + eV)}P(E-E' 
eRr Jo Jo 

By negating the integral variables in Eq. (|S-10|) and employing the symmetry of the Fermi function f(—x) = 1— f(x), 
of the BCS density of states ns(-x) = ns(x), and of the P(E) function P(—E + E') = P(E — £"), we obtain 

I{y) = ^lo dEn os(E)[fN(E-eV)-f N (E + eV)}, (S-ll) 
where we have introduced the effective DOS 

/oo 
dE' n s (E')P(E - E 1 ), (S-12) 
-oo 

We stress that to obtain Eq. (|S-11I) which corresponds exactly to current through an ideal NIS junction without 
coupling to environment but with a modified DOS, we only assumed that the superconductor is at low temperature and 
the environment is at high temperature. In this case, it is not possible to distinguish only from the IV characteristics 
whether the observed DOS arises, e.g., from some intrinsic property of the superconductor such as lifetime broadening, 
or from the electromagnetic environment described by P(E). As reported in the original manuscript, we can make 
this distinction by engineering the environment thus showing that it gave the dominating contribution to the effective 
DOS. 
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EQUIVALENCE OF THE P(E) AND DYNES DENSITY OF STATES 

The Dynes DOS 0, Q is given by 

E/A + i-y 



V(£/A + z 7 ) 2 -l 



(S-13) 



from which the BCS DOS ns is obtained in the limit 7 — > 0. Thus ns(E ± i-fA) = rig(E), where the sign can be 
chosen freely. To show that the effective DOS arising from the environment in Eq. (jS-81) is equivalent to the Dynes 
model with 7 = 2irRkBT cmj / (Rk A), we write the inverse Fourier transform of rig as 



/OO 
n s [u-i sign(t) 7 A]e^*dw/\/2^ = e — r^l«l jr-i _ 
-00 



On the other hand, 



= T [J- 1 (ng )] = J" [e-^F- 1 (n s )] = ? (e^l) * n s = n eff , 



(S-14) 



(S-15) 



where * denotes a convolution integral and we have used the identity T 1 (/ *g) = T 1 (/) T 1 (g). Thus the effective 
DOS in Eq. (|S-12[) with the environment given by Eq. (|S-8|) is equivalent to the Dynes model. 
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